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The amplitude of displacive excitation of coherent phonons (DECP) in YBa2Cu307 observed in 
time-resolved spectroscopy increases by a very large factor of 15 between 20 K (below Tr) and 100 
K (above Tr). Microscopic modelling of phonon excitation reveals that these amplitude changes of 
the Aig mode phonons are related to the existence of the Van-Hove singularity (VHS) below but 
very close to Fermi energy. 

PACS numbers: 78.47. +p, 74.25.Gz, 74.25.Jb, 71.15.Fv 
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Femtosecond pump-probe (PP) measurements in 
YJia2C'us07-s indicated that ultrafast optical excitation 
causes oscillations in the photo-induced reflectivity^'^. 
For YBa2CM307, the oscillation frequencies determined 
by fast Fourier transformation (FFT) match with two 
Raman active Aig modes, corresponding to vibrations 
of Cu(2) (150 cm-i) and Ba (120 cm-^Y . Among all 
experimental facts, the strong growth of the oscillation 
amplitude between 20 and 100 K is most striking. Be- 



tween 20 and 100 K, the intensities of the 150 



and 



the 120 



mode change by a factor of 2 and 15, 



spectively. Mazin et al.^ interpreted this experimental 
evidence by a phenomenological model'* with some lo- 
cal density approximation (LDA) parameters and within 
the BCS framework of an isotropic gap. The different 
change of relative intensities of the Cu(2) and the Ba 
mode below and above Tc was interpreted as a conse- 
quence of the mixing of these two modes. The mecha- 
nisms of DECP above and below Tc are phenomenolog- 
ically different. Above Tc, the driving force for DECP 
is the usual interband transition and thus mixing of the 
two modes contributes little to the observed DECP. Be- 
low Tc, however, the DECP is additionally excited by the 
movement of the ions to their normal equilibrium posi- 
tions caused by the destruction of the superconducting 
pairs after absorption of the pump laser. In this case, the 
mixing of these two modes becomes much more impor- 
tant and results in a remarkable reduction of the Cu(2) 
mode. Although this phenomenological model captures 
the essential physics in the PP process (e.g, this model 
assumes the driving force for the DECP is proportional 
to the density of excited electrons), it does not consider 
the real physics on the microscopic level. 

In this Letter, we analyze the given experimental facts 
on a microscopic level combining the microscopic theory 
for DECP proposed by Kuznetsov et al.'' with the eight- 
band Hamiltonian model developed by Andersen et af . 
Our results indicate that the mechanism of DECP be- 
low and above Tc is microscopically different as claimed 
previously^and that the strong enhancement of the am- 
plitude of DECP below Tc is additionally caused by the 
existence of the Van-Hove singularity (VHS) approxi- 
mately 20 meV below the Fermi energy. We also study 
the influence of the gap symmetry on the amplitude of 



DECP. According to Ref.[5], the amplitude of the DECP 
for an Aig mode should satisfy the following equation of 
motion: 



n 



= -2^^i.E^k<4ck>. (1) 
k 

Here Mt and < c^^^ci^ > are the electron-phonon (EP) 
coupling and excited electronic density matrix in k space, 
respectively. 

It is known that three bands cross the Fermi energy 
(Ei?), two for electrons in the Cu02 plane and one of the 
CuO linear chain for the YJia2C'us07. Many important 
physical properties are due to the energy band near E^?. 
The eight-band model® is an effective, nearest neighbour, 
tight-binding model Hamiltonian for a bilayer Cu02 and 
describes very well the two "plane-bands" crossing E^?. 
The Hamiltonian in Ref. [6] includes the four a orbitals: 
tpl(Cud^2_y2), tp?,(Cus), V'|(02p^), tpl(03py) and four 
TT orbitals: Vi(02p^), Vi(03p^), ^j'^iCud^^) , ^j'siCud^y). 
This reduced Hamiltonian is a direct result of the so- 
called down-folding procedure. It is interesting to note 
that this model can describe the form of the extended 
saddle-points around the Fermi energy in the Brillouin 
zone (BZ). The corresponding VHS in the density of 
states (DOS) is regarded as the reason for enhancement 
of the superconductivity transition temperature^. Next, 
we will show in Fig.l the calculated the EP coupling ma- 
trix in the frame of the above bilayer eight-band tight- 
binding(TB) model. A more general description of calcu- 
lating EP can be found in Refs.[8,9]. After the standard 
transformation*'^ into the momentum representation and 
diagonalization of the TB Hamiltonian of eight-bands in 
Ref. [6], we obtain the reduced Hamiltonian including 
the EP coupling and get the EP coupling matrix for the 
Ala mode®: 



M^,^(k,q= 0) = eE, 



1 
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:[ J2 l^l,iTM-i^b)/E. J2 I'l.iTM (2) 



s j'=3,4,5,6 



s j' = l,2,7,8 



1 



where //s j (7,k) are the coefficients of the diagonalizing 
the H^{s) Hamiltonian contributed by different atomic 
orbitais (j), bottom (s=0) and top (s=f) layer, respec- 
tively. Other notations are same as in Refs.[8,9]. 
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FIG 1. The evaluated EP coupling matrix M^j_^(k) 
for the Aig mode of the top layer (s=l) in the bilayer 
Cu02 system, as a function of k in the 2D lattice. The 
parameters for calculating Eq.(2) are taken from the di- 
agonalization of the Hamiltonian in Ref.[6]. 



relaxation of electrons and shows substantial changes of 
the PP signal for Y Ba2C'us07-s samples with different 
oxygen content but it has no effect on the excitation 
of the phonons in Cu02 plane. We will not consider 
this transition in this paper. Thus, only electrons in 
the vicinity of the S-point can be excited and the cor- 
responding integration conditions for the right side of 
Eq(l) would be: ( where rj can be choosen as ±7r/2a ) 
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This result is strongly supported by the density of 
states(DOS). Figure 3 shows the total (DOS) calculated 
from the bilayer eight-band TB model. The DOS ex- 
hibits two Van-Hove singularities (VHSs) below E^?: one 
very close to E^? ( approximately 20 meV below ) and an- 
other much farther away. In the eight-band TB model, 
the wave functions are divided into odd and even with 
respect to the mirror plane through the yttrium atom in 
the unit cell oiY Bu'^CusOt. The corresponding VHS for 
an even band is indicated by the solid line and the DOS 
for an odd band by the dashed line. Further calculations 
indicate that the main contribution to the DOS around 
the S point near the Fermi energy Ep comes from odd 
rather than from even bands. Therefore, only electrons 
described by the dashed line will thus contribute to the 
DECP in the normal state. 
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FIG 2. The IDA band structure of YBasCMsOr (from 
Ref.[6]). For 2 eV photoexcitation, the inter-plane-band 
transition from saddle-points around point X or Y are 
forbidden in the normal state because the transition en- 
ergy is much larger than 2 eV. 
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FIG 3. The total density of states (DOS) calculated 
from the eight-band TB model (solid line). The dashed 
line indicates the contribution of the even-bands ( see 
text ). 



The band structure of YBa2CM307 in the normal state 
obtained from LDA calculations is shown in Fig. 2. The 
saddle points corresponding to the VHS are predicted 
around the X and Y points as confirmed by photoe- 
mission results"'^"'^ . It can be seen, at normal state, that 
the optical inter-plane-band transitions (shown by the 
arrows in Fig. 2) around X or Y are impossible for 2eV- 
photoexcitation because the transition energy is much 
higher than 2 eV. Another direct optical transition would 
be possible from the saddle-point around Y to the chain 
band C. This transition plays an important role for the 



The interband transition ( i.e. electron/photon inter- 
action alone) hardly can explain^ the substantial change 
of the refiectivity below Tg. In the superconducting state, 
however, breaking of Cooper pairs which may be also ac- 
complished via electron/electron or electron/phonon in- 
teraction represents an alternative mechanism for DECP. 
Although direct interband photoexcitation of unpaired 
electrons and Cooper pairs at the saddle points around 
X or Y is still impossible, effective Cooper pair breaking 
can occur via collisions with unpaired carriers excited by 
the 2 eV photons at the S point. Additionally the fast 
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relaxation of the photo-excited paired- and unpaired- 
electrons around the S point in the BZ results in the 
generation of a large density of phonons whose energy 
is comparable or larger than the superconducting gap. 
Breaking of Cooper pairs at the X and Y point via elec- 
tron/electron or electron/phonon scattering can explain 
the drastic increase of the the driving force for DECP in 
the superconducting state because of the existence of the 
VHS at these points in the BZ. As shown in Fig. 3, the 
VHS around Ep comes from the even band. The excited 
electrons will be described by the solid line in Fig. 3 in 
the superconducting state. The limits for the integra- 
tion of the right side of Eq.(l) in the BZ are now given by 

(b): k e BZ 

Next we will try to estimate the DECP in the super- 
conducting state. For simplicity let us assume at the 
beginning that the gap is isotropic. Then the excitation 
energy of quasiparticles in the superconducting state is^^ 

^k,. = ^4,.+A(k)2 (3) 

The actual number of broken pairs in general would 
be less than the total number of superconducting pairs. 
However, the LDA estimation in Ref.[3] shows that the 
actual number of broken pairs will be very close the to- 
tal number of superconducting pairs in the ground state 
for the typical excitation intensities applied in PP ex- 
periments. Therefore we can take ^j^. ^ (5(^k,s — E) to 
calculate the number of excited electrons and their con- 
tribution to the driving force for DECP in the supercon- 
ducting state. In addition to pair breaking, the pump 
pulse also excites unpaired electrons whose contribution 
is assumed to be identical with the normal state. Thus 
according to Eq.(l) and the integration limits a and 6, 
the ratio of the amplitudes (ROA) of DECP below and 
above Tc will be 
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FIG 4. Numerically calculated ratio Da-^ s/Da-^ 
(see Eq.(4)) as a function of the Fermi energy shift 5 
normalized to the energy of Aig mode. For the d-wave 
case dj;2_y2 and d^y gap symmetry are considered. 
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The Van-Hove scenario will easily explain the occur- 
rence of a high superconducting transition temperature 
for several pairing mechanisms^'*. Indirect evidence for 
the existence of such a singularity has been reported in 
studies of the specific heat, thermopower, isotope effect, 
etc. Recently, direct evidence for a saddle-point, related 
to the Cu02 planes and located about 19 meV below the 
Fermi energy E^? at the Y or X point in the BZ, has 
been obtained in high energy resolution, angle-resolved 
photoemission spectroscopy (ARPES) on YBa2CM408^^ • 
It is worthwhile to note, however, that optical studies in 
which the rather high k-sensitivity should be assumed 
have not led to any signature for the VHS in the cuprates 
so far. In the following, we will demonstrate that the 
15 times increased amplitude of the Aig mode yields 
additional evidence for the existence of the VHS approx- 
imately 20 meV below the Fermi energy. In the normal 
state, as we have discussed before, the driving force for 
the DECP is the 2 eV interband transition optically al- 
lowed only around the S point in BZ as shown in Fig. 2. 
Thus, the total number of excited electrons will be very 
small as shown in Fig. 3. In the superconducting state, 
however, in addition to the interband transitions, effec- 
tive pair breaking will contribute to the driving force on 
the right side of Eq.(l). More important, the contribu- 
tion from pair breaking will become much larger if there 
exists a VHS below but close to Ep. Consequently, the 
huge number of electrons located around the VHS will 
be involved in the PP process in the superconducting 
state. The ratio of amplitudes (ROA) of DECP below 
and above Tc can be calculated as a function of the shift 
5 of the Fermi energy normalized to the phonon energy 
of the Aig mode. This can be done by replacing Ep in 
Eq.(4) by Ep -\- 5. Positive shift means that the VHS 
is farther away from Ep. Figure 4 shows the sensibil- 
ity of the VHS position on the ROA. It can be noted 
that the calculated results agree with PP measurement 
quite well in the case that the VHS is around the Fermi 
level. Consequently, it is impossible to observe an in- 
crease of the DECP amplitude below Tc by a factor of 
approximately 15 if the VHS does not exist or is far 
away from the Fermi energy. Although Dagotto et al 
proposed explanation of the fiat bands as coming from 
the correlation effects rather than band structure*®, the 
origin of the fiat bands is still unclear. In this Letter, we 
showed that the VHS described well by one-electron ap- 
proximate LMTO calculation can explain quantitatively 
the changes of DECP of Aig phonon in the PP measure- 
ments. Actually, as shown by Eq.(l), VHS scenario, no 
matter whether it comes from strong correlation effects 
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or band structure, can explain qualitatively the optical 
data in the time-resolved spectra. 
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different. The influence of the gap symmetry on the am- 
plitude of the Aig mode is rather small but becomes very 
large for Big or modes. The model calculations sug- 
gest that for Aig mode, R0A~15, while for Big mode, 
ROA is predicted as ~110 for dj;2_y2-wave and ~20 for 
s-wave symmetry. 
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FIG 5. Numerically calculated results of equation (4) 
as a function of the Fermi-energy shift S normalized to 
the phonon energy kujo of the Big for either s or d-wave 
( with dj;2_y2 and d^y gap symmetry) coupling. 

Finally the infiuence of the gap symmetry (GS) is now 
considered. This can be done by simply inserting the 
anisotropic gap equation A(k)=AoC(k) into Eqs.(3) and 
(4). GS is dx2-y2 with C(k) = (coskx — cosky)/2 and 
dxy when C(k) = sink^sinky. For d^y GS, the nodes 
will appear at the X and Y points (i.e. around the 
saddle-points of the VHS). In this case the carriers at 
the VHS form no Cooper pairs and thus the VHS cannot 
contribute to the driving force for the DECP. Therefore 
the result of the PP experiment excludes d^y GS. The 
numerically calculated curve shown in Fig. 4 indicates 
that the ROA of the Aig mode increases by a factor of 
16 for dj;2_y2 pairing, which is almost same as the value 
of 17 predicted for an isotropic gap. However, for the 
Big or B2g mode, the ROA of DECP is even much more 
sensitive to the symmetry of the gap. In Fig. 5, we show 
the same calculations for the Big mode. The ROA of this 
mode will increase by a factor of 110 for dj;2_j,2 pairing, 
as compared to a factor of 20 for the isotropic pairing. 
Again, ROA for d^y is also much smaller than value for 
other gap symmetry. Unfortunately, no experimental 
data have been reported for the DECP of the Big or B2g 
mode so far. We would like to emphasize, however, that 
according to our theoretical model studies of 5ig or B2g 
modes in optical experiments may provide a powerful 
tool to discern the gap symmetry. 

In summary, we have shown that the VHS scenario 
can lead to very large amplitude of the DECP for 
YBa2CM307 below Tg. The mechanisms for excitation 
of the DECP below and above Tc are microscopically 
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